Peierls Instability Induced Ferromagnetic Insulator at Orbital Order Transition 
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The origin of ferromagnetic insulating state of Lay/gSri^gMnOa is investigated. Based on the tight- 
binding model, it is shown that this state can be attributed to the Peierls instability arisen from 
the interplay of spin and orbital ordering. The importance of the hole-orbiton-phonon intercoupling 
in doped manganites is revealed. This picture explains well the recent experimental finding of the 
reentrance of ferromagnetic metal state at low temperature [Phys. Rev. Lett. 96, 097201 (2006)]. 
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I. INTRODUCTION 

Ri-zAzMnOs (R being rare-earth ions and A for di- 
valent ions, e.g. A= Ca, Sr, Ba or Pb) have been in- 
tensively studied for more than a decade because of its 
rich physics, such as unusual colossal magnetoresistancc 
(CMR) 1 , interesting ferromagnetic (FM) and antiferro- 
magnetic (AFM) phases, and charge and/or orbital or- 
dering. Despite intensive research and progress made so 
far in understanding the system, the origin of the fer- 
romagnetic insulating (FMI) phase observed in lightly 
doped L^/gSrjygMnOj 2 - is still controversial and un- 
known. CMR can be qualitatively understood in terms 
of the celebrated double exchange (DE) model 3 -, which 
means the simultaneous appearance of ferromagnetism 
and metallic behavior. Thus, FMI is not compatible with 
the DE model, at least not in a direct way. Although it is 
generally recognized that the FMI state may relate to the 
interplay of charge, spin and orbital degrees of freedom, 
no concrete picture is available. Finding such a picture 
is the focus of this paper. 

It shall be useful to first summarize the experimental 
observed phases of Lay/gSrx/gMnOa solid. Lowering the 
temperature, the first structural phase transition from or- 
thorhombic to monoclinic orders occurs at Tjt ~ 283K. 
The system is in the orthorhombic paramagnetic insulat- 
ing phase (phase I) above Tjt, and monoclinic param- 
agnetic insulating phase (phase II) below Tjt but above 
ferromagnetic (FM) transition temperature Tc ~ 183K. 
The second structural phase transition from monoclinic 
to triclinic orders occurs at Too ~ 150K, separating the 
monoclinic ferromagnetic metallic phased, (phase III) 
from the triclinic FMI state (phase IV) having a super- 
structure with unit cell 2a c x 46 c x 4c c (a c x b c x c c 
is the unit cell of the high-temperature orthorhombic 
phase)£. The FMI state of phase IV is supported by 
thermal activated exponential T-dependence 2 -^ of the re- 
sistivity that increases rapidly with the decrease of tem- 
perature for T < Too j an d the metallic nature of phase 
III is confirmed by the fact that the resistivity decreases 
with temperature below Tcr 2 -. The system undergoes an- 
other metal-insulator transition from FM insulator to FM 



metal (V) at temperature Tr ~ 30K. Phase V is a reen- 
trance of FM metallic (FMM) state through the FMI 
phase in Tr < T < Too? newly discovered by nuclear 
magnetic resonance technique 8 . The crystal structure be- 
low Tr is still unknown. Table 1 lists all five phases, char- 
acterized by four critical temperatures. So far, phases 
I-III have been widely studied and well understood in 
terms of crystal structure, magentic and transport prop- 
erties. 



TABLE I: Phase transitions of 
P — Par amagnet ic, FM— Ferromagnetic, 
M-Metallic 



La 7 / 8 Sr 1/g Mn03: 
I— Insulator and 



Phase 


Temperature 


Structure 


Magnetism 


Transport 


I 


T > Tjt 


Orthorhombic 


P 


I 


II 


T c < T < Tit 


Monoclinic 


P 


I 


III 


Too < T < T c 


Monoclinic 


FM 


M 


IV 


T R < T < Too 


Triclinic 


FM 


I 


V 


T < Tr 


Unknown 


FM 


M 



In order to explain the FMI state of Lay/gSr^gMnOs, 
many different models were proposed in the literature, in- 
cluding charge polarons order—, checkerboard-like charge 
order—, orbital order without charge order—, etc. Most 
of them are not consistent with the recent resonant X- 
ray scattering (RXS) experiments, which observed an or- 
bital polaron lattice (OPL) in phase IV—. This obser- 
vation is the basis of FMI— orbital polaron model, pro- 
posed by Kilian and Khaliullini 3 -. Orbital polarons can 
be viewed as charge carriers (holes) dressed by the inter- 
connected orbital states. At light doping magnanites, 
one Mn + site is surrounded by six neighboring Mn 3+ 
sites, equivalent to one hole (with e g orbits unoccupied) 
surrounded by six (occupied) e g orbits. The strong hole- 
orbital coupling will polarize those e g orbits and make 
them point towards the hole to minimize the interaction 
energy as well as the kinetic energy™ The formation 
of OPL in La7/gSr!/gMn03 emphases the importance of 
hole-orbital interaction for the FMI stated—. 

Obviously the genuine mechanism of the FMI state 
should explain following experimental features: 1) the 



2 



measured p—T curves 2 ^; 2) the giant phonon softening 113 ; 
and 3) the reentrance of the FMM state (IV-^V)£. It is 
difficult for orbital polaron model to explain feature 1) 
and 3). The reasons are as follows. Firstly, the orbital 
polaron cannot reproduce the measured p — T curves of 
Lai-zSrzMnCb at T < Too (cf. Fig.[3]). Polaron hop- 
ping model yields p — po exp(T a /fcsT), where E a is the 
activation energy of polarons, from which the binding 
energy (Eb = 2E a ) of a polaron can be deduced. The fit- 
ting to experimental data^ 2 - with hole-polaron hopping 
model gives E b ~ 0.06 — 0.12eV, about an order of mag- 
nitude smaller than the theoretical value (Eb ~ 0.6 eV in 
Refll3T). Secondly, the observed reentrance of the FMM 
state at lower temperature is inconsistent with polaron 
model since polaron hopping becomes more and more 
difficult when the temperature is lowered. Therefore, the 
orbital polaron itself is not sufficient to explain the FMI 
state and the origin of the metal-insulator transitions at 
Tqo remains intriguing and controversial. 

The above difficulties can be removed by including 
electron-phonon ( e-ph) interaction in the simple electron- 
orbit coupling model. We argue that the mechanism 
of the FMI state is the e-ph coupling induced Peierls 
instability 16 that opens an energy gap at T < Too- In 
this paper, we shall show that the Peierls instability can 
explain the observed reentrance of FMM state quantita- 
tively, and the other experimental features at the quali- 
tative level. 



II. THEORY 

Our theoretical model is based on the quasi-one- 
dimensional (ID) confinement of the motion of holes in 
Lar/gSrjygMnOs, which has been demonstrated by fol- 
lowing experiments: 1) the convergent-beam electron 
diffraction and selected-area electron diffraction seeing 
the superstructure (2a c x 4b c x 4c c ) of the FMI phase 6 
(anisotropic 3D motion); 2) the resonant X-ray scatter- 
ing showing an alternation of hole-rich and hole-poor 
planes in c direction-, which confines holes to move 
in the two-dimensional hole-rich a — b planes (2D mo- 
tion); and 3) the resonant X-ray scattering also reveal- 
ing the formation of OPL 12 , which further confines holes 
to move along the one-dimensional charge stripes in a 
direction (ID motion, see Fig.Q] and the following ex- 
planations). It is known that the quasi-lD confinement 
originates from the orbital order and its induced effec- 
tive coupling between adjacent Mn 3+ and Mn 4+ sites. 
The insert of Fig.Q] depicts two possible configurations 
of one Mn 4+ — 2 ~ — Mn 3+ unit with different e g orbital 
occupations on the Mn 3+ site. Configuration II corre- 
sponds to e g orbit pointing towards the hole along the 
axis of oxgen 2p orbit while configuration I corresponds 
to e g orbit pointing towards other directions. Due to or- 
bital anisotropy, configuration II results in a maximized 
overlapping of wavefunctions of occupied e g and 2p or- 
bits. That induces two direct effects: 1) lower energy 
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FIG. 1: (color online) The orbital order induced quasi-lD con- 
finement of the motion of holes in a direction. The dashed 
lines denote the weak interchain couplings in 6 and c direc- 
tions. The arrows near the O 2- ions indicate the displacement 
of the oxygen sublattice. The insert shows the orbital order 
induced effective coupling between adjacent Mn 3+ and Mn 4+ 
sites mediated by the oxygen 2p orbit. Configuration I and II 
correspond to different occupation of e g orbitals on the Mn 3+ 
site. A or j, is the splitting of energy between I and Ili 3 -. 



of configuration II than that of I with an energy dif- 
ference A or ir£; 2) effective Mn 3+ — Mn 4+ coupling with 
considerable transfer integration (t ^ 0). As a result, 
configuration II is stable and favors the motion of holes 
along a-axis. In the hole-rich planes shown in Fig.[TJ con- 
figuration II periodically repeats itself along a-direction 
thus forms the quasi-lD Mn 3+ -0 2 ~-Mn 4+ -0 2 ~ chain- 
like pathways of hole transport. Please be noted that 
weak interchain couplings (t±) in b and c directions also 
exist, resulting from the hole transfer between unoccu- 
pied and occupied e g orbits^ 2 -. Why the orbital order in- 
duced quasi- ID confinement can only be observed around 
x = 1/8? Qualitatively speaking, the orbital disorder- 
order transition takes place at light hole doping provided 
that the Jahn- Teller phonons and superexchange pro- 
cesses mediate an effective coupling between orbits on 
neighboring sitesi^. 

Now we are ready to present our model Hamilto- 
nian to describe the quasi-lD motion of particles (elec- 
trons/holes). As shown in Fig.Q] the particles mov- 
ing along Mn 3+ — 2 ~— Mn 4+ — 2 ~ chains will simul- 
taneously couple to the orbits on the manganese sites 
and lattice displacements (phonons) on the oxygen ones. 
Considering that the quasi-particles of orbital degree of 
freedom are bosonic orbiton o 1 i 18 , we can establish an 
electron-orbiton-phonon intercoupling model. For sim- 
plicity, only Mn sites are included in the model and the 
oxygen degrees of freedoms are integrated out, giving rise 
to the effective Mn-Mn coupling and the phonon modu- 
lation to it [cf Eq. (JTJ)] . Then, the tight-binding model 
can be written as (H = 1 and ks — 1 are assumed, and 
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they are restored in final results), 

H = t ^(ctcj+i + h.c.) + e?a+a g + ^ c j c i 2 J 
j q 3 

p j 

Qj = (Qj-i + Qj+t) 

Pj = (P s +i-P s ) (l) 

with 

q 3 = x;^^(o g +o± fl ) 

i>- = ^G p e^(6 p + 6± p ) (2) 
p 

where is the creation operator of a particle on site j 
(Rj is the atom position), et+ (6+) the creation opera- 
tors of a orbiton (phonon) with momentum q (p), e g and 
w p are the dispersion spectrum of orbiton and phonon 
respectively. The third term is the e-orb coupling be- 
tween a particle on site j and the orbitons on the nearest 
neighbor sites j ± 1, and the last term describes the e-ph 
coupling due to the phonon modulation of the particle 
hopping between sites j and j + 1— (see Fig.[T]). R q and 
G p are the e-orb and e-ph coupling constant respectively, 
which satisfy G*_ p — G v and R*_ q — R q to ensure the 
Hcrmitian character of Eq. |T|) . 

What are the different roles of the e-orb and e-ph 
interactions in Eq. (JXJ) ? One can see that the orbiton 
couples to the on-site charge density rij, which pro- 
motes the formation of Holstein-like small polarons in 
the strong-coupling limit 2 ^. That is the physical ori- 
gin of the formation of orbital polaron o 12 ' 13 . The essen- 
tial effect of Holstein-like small polarons is reducing the 
electron bandwidth (or enlarging the electron effective- 
mass). This interaction predicts a transition from band 
transport to thermal-activated hopping of polarons at 
high temperature, but it can not open an energy gap 
to induce the metal-insulator transition at low tempera- 
ture. However, the e-ph interaction in Eq. ([T]) can take 
the role of opening the energy gap, since the phonon 
modulation of the electron transfer predicts the Peierls 
instability in quasi-lD systems at low temperatur o 16 ' 21 . 
For Lay/gSi'i/gMnOs, when the orbital order effectively 
confines the particles to move along quasi-lD pathways 
(see Fig.[T]), the strong e-ph interaction distorts the oxy- 
gen sublattice to make O 2- ions displace towards Mn 4+ , 
which enhances the polarization of orbits and helps 
to stabilize the orbital order (see Fig.[T|). There thus 
forms the quasi-lD charge-density-wave (CDW) state 
in manganese sublattice accompanied with bond-order- 
wave (BOW) state in oxygen sublattice with twice the 
period of the original lattices (dimerization) . One anal- 
ogous material with similar CDW/BOW ground state is 
the halogen-bridged mixed- valence metal complex which 
has been well studied with this kind of e-ph interaction 
in the literature 2 ^. 



At T < Tq, the ferromagnetic t2 g spins of Mn ions 
exclude the opposite spin to occupied e g orbitals due to 
the Hund's rules, i.e. only carriers with the majority- 
spin contributes to the charge transport. Before coupling 
to orbitals and phonons, the Mn 4+ — Mn 3+ transfer [the 
first term in Eq. ([I])] in majority-spin subspace conduces 
a half-filled metallic energy band and the bare particle 
spectrum £® can be written as 

e k =i° k +s° k 0) 

with 

7° = — 2t cos k x a 

5® = — 2try(cos k y b + cos k z c) (4) 

where the effect of weak interchain coupling (t± = rjt, 
t] <C 1) has been included in r\ is a dimensionless 
parameter denoting relative strength of the interchain 
coupling to intrachain one. 

The orbitons and phonons then couple to the tight- 
binding particles, as described in Eq. ([1]), and change 
above metallic band to a insulator one. The simultaneous 
couplings between two bosons to one fermion are techni- 
cally difficult to solve, even though the direct orbiton- 
phonon coupling has been ignored in Eq. |T]). Thus we 
make two approximations for e-orb coupling in follow- 
ing calculations. The first one is the g-independence 
e q (= e) and R q {= r), based on the theoretical result 
that the orbiton is almost dispersionless at strong e-orb 
couplings 23 '. The second one is the perturbation theory 
of small polarons 1 ^ at strong e-orb couplings, which con- 
duces 

& - 4° exp(-r 2 / £ 2 ) - - = A$ - B (5) 
e 

where A denotes the electron bandwidth reduction (or 
effective-mass enhancement) and B the energy shift due 
to strong e-orb couplings. The above approximations 
have captured the essential physics of e-orb coupling: 
electron effective-mass enhancement (A) and energy shift 
(B), so that it should produce reasonable results in the 
limit of strong e-orb interactions. Please be noted that 
above approximations are only valid at finite e to avoid 
B unphysically large. After electron coupled to orbital 
degree of freedom, 

6c = 7fe + h (6) 

where 

Ik = A^ 

4 = A5° k -B (7) 

Let us make some further comments on the e-orb in- 
teractions concerned in the present work. Our model 
Hamiltonian describes electrons coupled to two bosonic 
fields (orbitons and phonons), which is formally simi- 
lar to the model of electrons coupled to two kinds of 
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phonons for cuprates^. The key feature of orbitons dis- 
tinguishing them from usual phonons is the strong spatial 
anisotropy due to the symmetry of e g electron wave func- 
tions, which has been qualitatively considered in building 
up our quasi-lD Hamiltonian [see Fig.[T]and Eq. (JT}]. The 
quantitative description of orbitons by means of pseu- 
dospin operators, which may be required for more gen- 
eral studyi^, has been beyond the scope of this paper. 
The parameter characterizing the orbital ordering come 
into play due to orbitons is A = r/t. In our calcula- 
tions, A Q ~ 0.2 — 0.65, which is realistic to manganites 
and consistent with literatures^. 

Particles with new energy spectrum in Eq. ([6]) then 
couple to phonon degree of freedom. Now, we can an- 
alytically solve the Hamiltonian by transforming it into 
momentum space, 

H = H ^4 c k + w p 5 p b P + H 9pkct +p c k (bp + b± p ) 

k p k } p 

where g pk is the e-ph coupling constant in momentum 
space and its dependence on k is to be neglected in what 
follows (g p — g*_ p oc M p ), as in the literature^. That 
Hamiltonian for ID systems indicates a Peierls super- 
structure which, described by introduction of the follow- 
ing anomalous average^, breaks translational symmetry 
of initial lattice: 

A = g 2kF < b 2kF + bt 2kF >+ (9) 

where angular brackets denote thermodynamic average 
that can be obtained by performing Gibbs average to 
Matsubara equations of motion for operators bx and b* K 
(K = 2Uf for Peierls phase transition)^, 



V(x) — Ae tKx + A*e~ tKx . The anomalous Green's func- 
tion F can be derived from the Gorkov equations for 
Matsubara Green's functions (limit to first order in V) 
under the "nesting condition" : jk-K — 7fc-2fe F = 



-7A 



G(k,e n ) = G Q (k,e n ) + G (k,e n )AF(k,e n ) 
F(k,e n ) = G Q (k-K,e n )A*G(k,s n ), (13) 

which gives the following solutions, 

ie n - 5k + 7fc 



G(k,e n ) — — 



By replacing ie n with e, the new energy spectrum is de- 
termined by the zero of denominators (pole) of Eq. (|14p : 



e = 4 ± \Hi + A 



(15) 



Inserting Eq. ((7j) into Eq. (fl~5j) gives the energy gap E g of 
the system, 



E g =2A- 8Ar)t. 



(16) 



By inserting Eq. (fT4|) into Eq. (|T2|) and then perform- 
ing standard calculations, we obtain the self-consistent 
equation of A (set to be real), 



sinh v^p! 



2 V cosh cosh ^ V7 fc 2 + A 2 



(17) 



where g = g K y/2/cu K - 



(-^--uj K )<b K (r)> = £>(fr,T = -0) 

k 

{-^+"K)<b + _ K {T)> = -^F(fc,r = -0)(10) 



where F(kt) = — i < Tak(t)at_ K (0) > is the anoma- 
lous Green's function describing the elementary Umklapp 
scattering process k — K — > k. 

After Fourier transformation over Matsubara "time" 
to Eq. (fit)]) , we have: 



< b K + bt K > Um = - 



2g K uj K T 



E^- 6 ") ( n ) 



k . 7 1 



The condition for the Peierls phase transition is u> m = 
at K = 2kp, i.e., 



A = gx < b K + bt K >uj m =o^ 



2g^T 



Y,F(k,s n ) (12) 



In coordinate representation Eq. © describes Peierls de- 
formation potential characterized by the wave vector K: 



III. DISCUSSION 

So far we have derived the insulating state at T < Tq 
resulting from the Peierls instability at the orbital or- 
der transition temperature Too- With this mechanism, 
we first quantitatively explain the observed reentrance 
of FMM state. The key point is that a T-dependent 
interchain coupling rj(T) may induce two critical tem- 
peratures T c i and T c2 of the metal-insulator transition. 
The former corresponds to the ordinary Peierls instabil- 
ity and the latter to the reentrance of the metallic state. 
It is called stepped Peierls transition theory, with which 
Zhou and Gong explained the anomalous transport prop- 
erty in NbSeg^l. With this theory, we can elucidate the 
reentrance of FMM state in Lar/gSrjygMnOs. For strong 
e-orb coupling, the self-consistent equation of Tmi [de- 
rived from Eg (Tmi) — 0] can be expressed as 



A(Tmi) = 4Ar)(T MI )t. 



(18) 



In principle, rj(T) should increase with the decrease of 
temperature, in analogy with the role of pressure. Since 
the exact expression of rj(T) is hard to determined, we 
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take the leading linear term of its series expansion at low 
temperature (T < 7a), 

fj(T) = %(1 - aT/T A ), (19) 

where Ta denotes the critical temperature of the charge 
order [A(Ta) = 0], which coincides with T c i at 77 — [cf. 
Eq. (|16[) ] but is a little higher than T c \ at 77 > 0. 770 is 
the T-indcpcndent part of 77 and a a phenomenological 
parameter. In our calculations, we fix Ta = 180K by con- 
sidering the experimental T c i = Too ~ 150K. Please be 
noted that at least two parameters in the self-consistent 
equation [Eq. (IT71) ] dependents on each other after Ta is 
fixed. 

We then self-consistently solve Eq. (flTt - pl?)) and il- 
lustrate the change of T c i and T C 2 with 770 in Fig. [2] 
The parameters are chosen as r = 0.2t, g = l.ly/t, 
e = 0.6t, a = 0.5 and t = 0.4 eV. As shown in Fig.H 
our calculations can yield T c i ~ 150K and T C 2 ~ 30K 
at 770 ~ 0.028, comparing well with the experimental ob- 
servation in Ref. HI as indicated by the big arrows in the 
figure. No reentrance of metallic state occurs at small in- 
terchain coupling (770 < 0.024), characterized by T c2 = 0. 
T c i decreases slightly with increasing of 770 in this limit. 
Nonzero T C 2 emerges at about 770 ~ 0.024. Further in- 
creasing interchain coupling will induce rapid increas- 
ing T C 2 and rapid decreasing T c i meanwhile. When T c i 
and T C 2 coincide into one point at large interchain cou- 
pling, the Peierls phase transition disappears. This is just 
the intcrchain-coupling-induced derealization of quasi- 
1D states, as occurs in conjugated polymers^. Since the 
value of t for e g electron system is believe to be 0.3 ~ 0.6 
e V 13 i 29 , our theory quantitatively explains the reentrance 
of the FMM state with consistent parameters in the lit- 
erature. 

In a recent experiment on the (Oll)-oriented 
Lar/gSriygMnOs films, Chen et al. measured an energy 
gap opening at T ~ 190K and increasing to ~ 0.16 eV 
at T = 120K 30 . The absent physics mechanism in their 
experiment has been already elucidated in this paper, 
that is, the observed energy gap results from the Peierls 
instability at orbital order transition. To further clarify 
this point, we calculate the temperature dependence of 
the energy gap self-consistently, and show the result in 
the insert of Fig. [5] together with the experimental data. 
Some of above parameters are adjusted to fit Chen et 
aVs experiment: T A = 200K, r = 0.64* and 770 = 0.02. 
As shown in the figure, our calculated E g — T curve is in 
good agreement with the experimental data. 

We then systematically study the effect of the e-orb 
coupling on the Peierls phase transition, and find that 
large r increases T c \ and decreases, in meanwhile, T C 2. 
As a consequence, strong e-orb coupling will enlarge the 
phase area of the FMI state in Fig.[5J Sufficient strong 
e-orb coupling will annihilate the reentrance of the FMM 
state by decreasing T C 2 to zero (the figure not shown). 
That confirms the orbital order being in favor of the 
Peierls insulating phase. Please be noted that the Peierls 
instability can not take place without the orbital order 







T 1 1 






■ 


FMT 








L f 1V11V1 


0.15 




a Exp. 












— Theo. 






j \ 


r> 0.10 












ta 












"To.05 






















i h 












./ CI 


0.00 












120 


140 160 180 


200 








r(K) 












J 1 1 






'FMM 



0.01 0.02 0.03 



FIG. 2: (color online) Two critical temperatures of metal- 
insulator transitions (T c i and T C 2) as a function of interchain 
coupling strength 770. The large dots are the coincident points 
of Tci and T C 2- The reentrance of the ferromagnetic metal 
state is shown by the big arrow. The insert shows the temper- 
ature dependence of the energy gap. The line is the theoretical 
result while the dots are experimental data from Ref. [3(J See 
text for details of the parameters. 
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FIG. 3: (color online) The fits of the experimental measure- 
ments of p by p = po exp(To/T) 1//2 . The experimental data is 
from Ref. Qj (x = 0.125) and Ref. Q (for others). 



induced quasi-lD confinement in Lay/sSriygMnOa. 

With our theory, we can also explain the other exper- 
imental features listed in the introduction at the qual- 
itative level. The Peierls instability opens an energy 
gap [see Eq. ([I6"]l] and makes Lay/gSri/gMnOa undergo a 
metal-to-insulator transition at Too- As a consequence, 
the charge carriers transport only through the localized 
gap states at T < Too, which will trigger the variable 
range hopping (VRH) 31 . That picture is confirmed by 
the the good fit of experim ental p—T curvei^ with VRH 
mechanism: p — po exp( 1k /T /T) in Fig. [3] The exponent 
1/2 results from a parabola density of states (DOS) with 
zero-DOS at the Fermi level Ep (called Coulomb gap). 
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As regards the the giant phonon softening of the Mn-0 
breathing mod&i^, we comment that it results from the 
celebrated giant Kohn anomaly at Peierls instability, i.e. 
the suppression of phonon frequency at p ~ 2fci^2. 

IV. SUMMARY 

In summary, we have demonstrated that the FMI state 
of Lai_ x Sr x Mn03 originates from the electron-phonon 
coupling induced Peierls instability when the orbital or- 
der confines holes to move along quasi-lD pathways. 
With this picture, the reentrance of the FMM state has 
been well explained quantitatively. The other experimen- 
tal features of the FMI sate, such as the temperature de- 



pendents of resistivity and the giant phonon softening, 
have been also understood at the qualitative level. Our 
theory supports the belief that the intercoupling of hole- 
orbital-phonon is critical in understanding the electronic 
properties of doped magnanites. 
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